Stable Spatial Langmuir Solitons. 
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We study localized two- and three-dimensional Langmuir solitons in the framework of model based 
on generalized nonlinear Schrodinger equation that accounts for local and nonlocal contributions 
to electron-electron nonlinearity. General properties of solitons are investigated analytically and 
numerically. Evolution of three-dimensional localized wave packets has been simulated numerically. 
The additional nonlinearities are shown to be able to stabilize both azimuthally symmetric two- 
dimensional and spherically symmetric three-dimensional Langmuir solitons. 
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Starting from the early 70-th, the fundamental role of 
Langmuir solitons (LS) in strong plasma turbulence is 
commonly known |2| • Formation of LS are connected 
with the plasma extrusion from the regions of strong 
high-frequency electric field and trapping of plasmons 
into the formed density well (caviton). However, in the 
previous simplified theoretical models, two-dimensional 
(2D) and 3D solitons occurs to be unstable with respect 
to a collapse: above some threshold power, the size of 
the structure shrinks infinitely, forming density singular- 
ity at finite time ilj. Physically, when the size of col- 
lapson approaches few Debye radii r^i, it should damp 
rapidly which may result in fast plasma heating. Never- 
theless, experimental observations 0, 0,11, IE of quasi- 
2D and 3D Langmuir collapse demonstrate saturation of 
wave-packet's spatial scale at some minimum value be- 
ing of order (10 ^ iO)rD- It have been observed in Refs. 
lET^ that at times t > 50ujpl, where mpi is the ion 
plasma frequency, Langmuir wave packets show consid- 
erably slow dynamics (subsonic regime). To our best 
knowledge, these observations do not meet an appropri- 
ate theoretical explanation yet. 
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Various additional linear and nonlinear effects, such 
as higher-order dis pers ion |^ 9, 10], the saturation 
|nl 1 121 . nonlocal wave interaction 
llq. h 9[ ■ may arrest wave collapse 
both in 2D and in 3D |3, M> u2> IH Ql • As it is shown in 
|20l |. the local part of electron-electron nonlinearity (re- 
sulting from the interaction with the second harmonic) 
counteracts the contraction of wave packet. At the same 
time, the nonlocal contribution of the additional nonlin- 
ear term was omitted in [20] . though it is of great im- 
portance for sufficiently narrow and intense wave pack- 
ets. In this Letter we take into consideration both these 
extra nonlinear effects. As it will be shown below, the 
role of nonlocal nonlinearity is quantitatively even more 
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significant. The nonlocal nonlinearity is of great impor- 
tance not only when describing soliton formation in plas- 
mas 13, 21j 22], but also in the theory of Bosc-Einstein 
condensates or matter waves |3 0, El Q], and in 
the construction of an adequate continuum model of the 
electron- pho non interaction in discrete 2D and 3D lat- 
tices [1113 . 

The evolution of radial component of electric field 
strength i5 of a Langmuir 2D and 3D wave packet is 
described by the set of equations: 



.dE 



D 



d_ 

dr 



E\E\ 



dr 

3 £;ArlS|" 
2 487rmno'^p 



2 ?io 



0, 




167rM ' 



where azimuthal or spherical symmetry is supposed, r - 
radial coordinate, = ^ + ^T^^, '''d - Debye radius. 
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electron plasma frequency, d is the number of space 
dimensions, m and M are the electron and ion masses, 
no, Tf, - background electron density and temperature 
respectively and cs = ^JT^/M is the ion sound speed, ni 
is the electron density perturbation. This equation set is 
valid if W/nT < (kro)''^, where W = E"^ /%tt, k being 
the effective wave number of the packet. In [20], the sim- 
ilar set of equations for Langmuir wave field was derived, 
however, the nonlocal nonlinear terms were omitted. 

We consider subsonic motions and neglect the terms 
with time derivatives in the second equation of above 
set. As result, this set is reduced to the single partial 
differential equation: 
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where the coefficients D, B, C, T are given by the ex- 
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pressions: 



D = -uj„rl,, B = 2 ^, 

2 ^ 32nMnocl 

Equation Q is closely related to nonlinear Schrodinger 
equation with modified linear term (proportional to D). 
The nonlinear part of this equation includes common cu- 
bic nonlinearity (term proportional to B) as well as non- 
local (term proportional to C) and local (term with F) 
parts of electron-electron nonlinearity. 

Equation ^ conserves the following integrals: the 
plasmon number 




and Hamiltonian: 




Other integrals (momentum and angular momentum) are 
equal to zero in the case under consideration. 
Let us show that the effective width reft 

r^^^l J r'\E\'dr, (4) 

of any wave packet governed by the Eq. is bounded 
from below in the most interesting case of self-trapped 
wave packets having negative Hamiltonian. If both C and 
r in Eq. are equal to zero, and H < 0, collapse occurs. 
However, any of additional nonlinear terms (proportional 
to C or r) prevents collapse in 2D as well as in 3D cases, 
or, in the other words, rcfr is bounded from below. 
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FIG. 2: Three-dimensional Langmuir solitons. (a) Plasmon number vs nonlinear frequency shift A for a = 1.5; (b) Effective 
soliton width vs plasmon number. The inset: soliton profiles corresponding to the points marked at the N{X) and Rcfi{N) 
curves. 



Let us define 



Eq. is reduced to 
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where all integrals are supposed to be finite. It is easy to 
sliow that 
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Since Hamiltonian Q is assumed to be negative, it is 
necessary that 



B I 1^1* dv>C I {V \E\' ] dr + T I r-^ \Ef dr. (6) 



Then we use the "uncertanity relation" of the form 

rlJ{w\E\^ydr>a{d)J\Efdr, (7) 

where a = 3/4 in 3D case and a = 1/2 in 2D case. One 
can see from inequality © that r| > T/B, thus, using 
Eq. (O one gets rf > T/B. Taking into account that for 
a localized wave packet r^g > r\, we finally obtain 
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Thus, if C 7^ (or if F 7^ 0), the wave packet can not be 
contracted to the size smaller than yJC /B (or ^JT /B). 
By introducing the dimensionless variables: 
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where a = 7/2. (In the limit a = equation © re- 
duces to equation (18) of Ref. ^^.) We will keep a as 
free parameter to study the impact of nonlocality on the 
properties of 2D and 3D Langmuir solitons and bearing 
in mind the other possible applications of model equation 
®. 

Soliton solutions of Eq. ^ have a form ^'(i?,T) = 
iJjiR) exp (iAr), where A is the nonlinear frequency shift 
of the soliton. The radial soliton profiles iIj{R) are found 
from the equation 
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We will start our consideration with 2D Langmuir soli- 
tons. Note the formal analogy between Eq. HlOfl for ra- 
dial electric field component in 2D case and NLSE for 
vortex solitons ^(R) exp (inup) with topological charge 
m = 1 (see, e.g. |lOj'). However, in spite of zero value of 
wave intensity at the soliton center, it has no phase dis- 
location and corresponds to the ground state with min- 
imum energy and zero angular momentum. In the con- 
sidered case of Langmuir wave structures, the phase does 
not depend on the radial coordinate, thus in Eq. IjlUI) 
ip{R) may be considered as real function. 
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Stationary states of Eq. were investigated ana- 
lytically and numerically. Analytical approach employs 
approximate variational method (see, e.g. [23) with the 
normalized trial function of the form 

E{r)^./N^ ^e^v(-ei'^), C - M^, (11) 

where the variational parameter /i characterizes the in- 
verse soliton's width: /i = 1/(V2 r-cff). We have returned 
to the nonscaled variables for variational analysis. The 
trial function Hll|l has a correct asymptotic near the soli- 
ton center. The Gaussian profile gives a good approxi- 
mation to soliton solutions, since, as it was argued in [l^ 
it represents an exact solution in the limit case of strong 
nonlocality. The variational parameter /ig corresponding 
to stationary solutions of the form Hll|l. is readily found 
after standard procedure (23j: 

^il = ^il{N- No)/N, (12) 

where y.1 = S/(4C -f 2r), iVo = IQt:D/B. Thus, 2D 
Langmuir solitons are formed only when some threshold 
value of plasmon number is exceeded (i.e. if > Nq). 
One can see that fJ.Q ^ a,t N ^ Nq, so t hat the effec- 
tive width is bounded from below: rcs > \/ (2C + T)/B, 
which agrees with previous general estimate given by Eq. 

The boundary- value problem described by the Eq. H1U|) 
with zero boundary conditions at the soliton's center 
{R — 0) and at infinity {R oo) was solved numerically 
by the shooting method. Solitons form two-parameter 
family with parameters a and A. For each given non- 
locality parameter a we present the dependence A^(A) 
known as "energy dispersion diagram" (EDD), which 
are given in Fig. ^ (a). Note that all solitons are sta- 
ble, which is similar to Vakhitov-Kolokolov criterion [l^ . 
since dN/dX > 0. To excite stable two-dimensional 
Langmuir soliton, the threshold value iVo of input power 
(critical plasmon number) should be exceeded. Vari- 
ational approach predicts normalized critical plasmon 
number to be Nq — IGtt. This is in a very good agreement 
with our numerical results, where Nq ~ 48.3. Moreover, 
this simple variational analysis gives a good description 
for all EDDs, and variational dependence N{X) fits bet- 
ter the computed one when the parameter a grows. It 
is illustrated in the Fig. ^ (a) for different a. One can 
see that at cr = 0, the shape of soliton differs sufficiently 
from the simple profile of the form Hll|l [compare Fig. 
(c) with Fig. n(d),(e)]. The typical soliton profiles are 
plotted in Fig. □(c)-(e). At the same plasmon number, 
the soliton width is larger for solutions of Eq. H1U|) with 
nonlocal nonlinearity (cr > 0) than for solutions with 
CT = 0. 

As it was stressed above, the effective width of Lang- 
muir wave packet in our model is bounded from be- 
low. Figure ^ (b) represents the effective soliton width 
as function of plasmon number. The Rcs decays mono- 
tonically and saturates at some nonzero minimum value 
^min- This minimum width increases when parameter 



a increases, as it was estimated above [see Eq. (jSJ]. In 
the dimensional variables, the minimum diameters of 2D 
Langmuir structures are of order lOro which is in agood 
agreement with those observed experimentally in |5(. 

We have performed similar analytical and numerical 
study for 3D radially-symmetric solitons. The corre- 
sponding A^(A) and Rcf[{N) dependencies are shown in 
Fig. [11(a), (b). In the 3D case, plasmon number depends 
on nonlinear frequency shift A non-monotonically, and 
the only stable soliton branch corresponds to dN/dX > 0, 
as it follows from Vakhitov-Kolokolov criterion. There- 
fore, in the 3D case, the threshold plasmon number Ncr 
needed to excite stable Langmuir soliton should be ob- 
tained from condition dN/dX = 0. Similarly to the 
2D case, the effective width Rcs of stable solitons de- 
cays monotonically when plasmon number N grows, and 
i?off — > i?min at N :s> Ncr- Wc havc found the threshold 
value to be of order Ncr ~ 850. 

We have investigated evolution of radially symmet- 
ric Langmuir structures numerically in the framework of 
Eq. (O. The integration time step was splitted into 
the linear and nonlinear parts, and both of them were 
performed in the radial coordinate domain. During the 
simulations, the conservation of integrals (|2J) and © has 
been verified, and if the change of any integral exceeded 
1%, the modelhng was stopped. We used two types of 
initial conditions: (i) perturbed stationary soliton solu- 
tion of Eq. H1U|) which was found numerically, and (ii) 
the Gaussian-like profile. At the boundaries we assumed 
that d [*(i?, t) /R] /dR = at i? ^ and *(i?, i) = at 
R — !■ oo. 

Stability properties of Langmuir solitons were found 
to be in a very good agreement with our analytical pre- 
dictions. Any localized wave packet having the number 
of quanta below the threshold value N < Ncr always 
spreads out, while an initially more intense packet with 
N > Ncr may form a localized structure after irradiat- 
ing a portion of plasmons. The dynamics of intense wave 
packet with N > Ncr is very intriguing. Figure |2| rep- 
resents the different evolution scenarios of perturbed 3D 
Langmuir solitons. The amplitude and effective width of 
slightly perturbed soliton solution belonging to the stable 
branch (with dN/dX > 0) oscillate in time as it is seen 
in Fig. O (a). However, when perturbed, an unstable 
stationary solution (having dN/dX < 0) manifests differ- 
ent evolution pattern: it can either monotonically spread 
out [see Fig. |3| (c)] or develop quasiperiodical motion 
[see Fig. |31(b)], depending on the initial perturbation. 
The latter quasiperiodical behavior resembles the oscil- 
lations between two opposite extreme states discussed in 
[l^lTgl . The initial increase of wave packet's intensity is 
followed by quasiperiodical oscillations with rather large 
amplitude. Thus, the additional nonlinearities actually 
prevent catastrophic collapse of any 3D Langmuir wave 
packets. The similar regime is known as "frustrated col- 
lapse" ^3 in theory of Bose-Einstein condensates with 
nonlocal nonlinear interaction. 

Our considerations are in a qualitative agreement with 
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FIG. 3: Typical examples of evolution of 3D Langmuir wave packets {a = 1.5). Left panels represent distributions of wave 
intensity |'I'(_R, Right panels show wave packet's amplitude ^'max vs time, (a) perturbed stable soliton with A = 0.1; 
perturbed soliton from unstable branch at A = 0.03 and different perturbations: (b) demonstrates quasiperiodical dynamics; 
(c) spreads out. 



experimental observations of Langmtiir wave structures 
in unmagnetized laboratory plasmas j^, Q ■ On the slow 
(subsonic) stage of evolution, which was observed at time 
t > 50ujpl, the field envelope almost stopped its ra- 
dial contraction and the size remained approximately un- 
changed (r is of order of few tens of td). It is important 
to note that according to |^ Q at this stage of Langmuir 
structure evolution, the beam-wave resonance is detuned 
and the beam decouples from the wave. Nevertheless, 
the common Zakharov set of equations describes field 
evolution only for rather short time {t < bOujpl) and 
then its prediction drastically deviates from experimen- 
tal facts because the wave packet becomes too intense 



(at peak E'^^^/AimT ^1) and too narrow. As it was 
demonstrated above, the electron-electron nonlinearity 
plays a crucial role in Langmuir structure behavior and 
it gives the qualitative explanation of saturation of the 
wave packet contraction. The dissipative effects such as 
Landau and transit-time damping seems to be negligible 
for structures of several tens Debye's radius characteris- 
tic size. Certainly, the wave packet may interact with 
high-energy electrons, therefore, after a long time (of or- 
der of several hundreds of Wp] ) , it may eventually damp 
due to wave absorption. Our considerations are valid at 
stage when nonlocal nonlinearities come into play and 
considerably slow down the wave packet contraction but 
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the wave absorption is still not essential. 

In conclusion, we have performed analytical and nu- 
merical studies of spatial (2D and 3D) Langmuir solitons 
in the framework of model based on generalized nonlin- 
ear Schodinger equation including both local and non- 
local electron-electron nonlinear ities. Their influence on 
intense and narrow Langmuir wave packets are of the 
same order, and both nonlinearities should be taken into 
account simultaneously. Any of them is able to arrest the 
Langmuir collapse. Both nonlinearities lead to the sat- 



uration of soliton width with an increase of the energy, 
but quantitatively the effect of nonlocal nonlinearity is 
more significant. All 2D Langmuir solitary structures 
are stable, while in 3D case, two soliton branches coexist, 
one is stable and the other is unstable. When perturbed, 
stable solitons demonstrate centrosymmetric oscillations. 
As for 3D solitons from the unstable branch, they may 
either spread out or oscillate quasiperiodically depending 
on perturbation applied. 
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